Electronic properties of the ID Frenkel-Kontorova model 
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The energy spectra and quantum diffusion of an electron in a ID incommensurate Frenkel- 
Kontorova (FK) model are studied numerically. We found that the spectral and dynamical properties 
of electron display quite different behaviors in invariance circle regime and in Cantorus regime. In 
the former case, it is similar to that of the Harper model, whereas in the latter case, it is similar to 
that of the Fibonacci model. The relationship between spectral and transport properties is discussed. 
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The Frenkel-Kontorova model describes a one dimen- 
sional chain of atoms/particles with harmonic nearest 
neighbor interaction placed in a periodic potential. It 
is a widely used model in condensed matter physics and 
nonlinear dynamics Qj. For instance, it has been used 
to model crystal dislocations 0] , epitaxial monolayers on 
the crystal surface || , ionic conductors and glassy mate- 
rials B, electron in a quasi ID metal below the Peierls 
transition S, charge density waves [||, Josephson junc- 
tions chains 0, and dry friction |sj . More recently, this 
model has been employed to study transport properties 
of vortices in easy flow channels and strain-mediated 
interaction of vacancy lines in a pseudomorphic adsor- 
bate system pJpf . Due to the competition between the 
two length scales, the spring length and the period of 
the on-site potential, the FK model exhibits a wealth 
of interesting and complex phenomena (see Jll]] and the 
reference therein). 

One of the most striking features of the FK model is the 
so-called transition by breaking of analyticity. It is shown 
by Aubry [Q that there exist two different ground state 
configurations for an incommensurate chain. The transi- 
tion from one configuration to another occurs when one 
changes the coupling constant K (see Eq.(||)). These two 
incommensurate configurations correspond to invariance 
circle and Cantorus of the standard map fl3|| , respec- 
tively. This transition is still discernible in a quantum 
FK model @. 

Although extensive studies have been done since its 
introduction, the FK model continues to attract active 
interests from different fields. Recent studies are concen- 
trated on phonon modes, because they are responsible 
for the heat conduction along the chain jn|[l(|. It is 
found that the on-site potential breaks the conservation 
of momentum, and makes the heat conduction in ID FK 
model obey the Fourier law Jl6| . 

However, the electronic property of the ID FK model 
is still unknown up to now. This topic is important 
from both fundamental and application points of view, 



because incommensurate and quasiperiodic structures 
appear in many physical systems such as quasicrys- 
tals, two-dimensional electron systems, magnetic super- 
lattices, charge-density waves, organic conductors and 
various atomic monolayers absorbed on crystalline sub- 
strates. As mentioned before, the FK model has been 
very successful applied in these systems. It is the pur- 
pose of this Letter to study this topic. 

The electron in a ID FK chain obeys the equation 



t (lp n+ i + lpn-l) + V n lp n = El\) n , 



(1) 



where to is a nearest-neighbor hopping integral which is 
set to 1 in this Letter, tp n is the amplitude of wave func- 
tion at the nth site and E the eigenenergy of electron. 
The on-site potential, V n — Xcos(x„), is controlled by 
parameters A and x^. {x„} is the configuration of an 
incommensurate ground state of the FK model, namely, 
{x^} minimizes the functional 
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(x n+ i - x n - a) 2 + K(l - cos(x n )), (2) 



where K is a coupling constant, and a the equilibrium 
distance between consecutive atoms, {x^ } is determined 
by an adjustable coupling constant K . 

In contrast to the Harper model and the Fibonacci 
model that have been often used to study electron in in- 
commensurate systems (see review articles |l7],[l8|], Ref. 
p9| , and the references therein), the FK model has two 
control parameters, A and K, and is thus more general. 
As we shall see below that the two control parameters 
not only make the problem more complicated but also 
enrich the physics of this model. Most electron proper- 
ties shown in the Harper model and the Fibonacci model 
can be recovered in the FK model. 

Any change in A and K will alter on-site potential V n , 
thus change the electron properties of the system. As 
it is well known that for each irrational number a/27r 
there exists a critical value K c separating the two con- 
figurations of ground state. K c — 0.9716354 • • • corre- 
sponds to the most irrational number, golden mean value 



1 



a/2n = (y/5— l)/2. In this Letter, we restrict ourselves 
to this particular value in numerical calculations as it is 
the mostly used one in the community. 

In order to study electron energy spectra, we first 
obtain ground state configuration for N atoms by gra- 
dient method Q for fixed boundaries, i.e., xq = 
and xn — 2ttN. a/2ir — (\/5 — l)/2 is approximated 
by a convergent series of truncated fraction: F n /F n +i 

(n = 1,2 ), where {F n } is a Fibonacci sequence. 

The number of atoms is chosen as N — F n . The electron 
eigenenergies are obtained numerically by the transfer 
matrix method. The transfer matrix Tn(E) is 

N 

T N (E)=l[T(n,E), (3) 

n=0 

where 



model does not have self-duality and the spectrum is 
asymmetric. It is known Jl7| that for a non-self-duality 
system, there exits a critical parameter above that all 
eigenstates are localized. Below this critical value, ex- 
tended, critical, localized states, and mobility edges co- 
exist. 

To study the eigenstates quantitatively, the Thouless 
exponent (^f{Ei)) and participation ratios (PR) are cal- 
culated for every eigenstate. The Thouless exponent is 
given by 

jiEi) = J \n\Ek - EMfydEj = ^ X>I S ' " E ^ 

(5) 

where p(E) is the density of states, and the participation 
ratio is, 




FIG. 1. Energy spectra versus A for electron in the 
chains with different values of K. (a), K = 0.4 < K c ; 
K = 1.6 > Kc- The chain length is N = 377. 

The allowed energies of electron satisfy condil 
\TrT N {E)\ < 2. Fig. illustrates the spectra as ft 
tions of A for different values of K. From Fig. |l|, 
can see that for K < K c , the energy spectra (see 
|l|a) is similar to that of the Harper model, namely 
band splits into subbands as A is increased from 0. A 
becomes larger than a critical value, energy levels t 
to repel each other. There are however differences fi uiii 
the Harper model. For instance, the Harper model has 
a good symmetry, such as self-duality, and the spectrum 
is symmetric about E — 0. All eigenstates are extended 
when A < A c (= 2), and localized when A > A c . Our 
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FIG. 2. The Thouless exponent y(E) for different values of 
K. (a), K = 0.4, curves 1, 2 and 3 correspond to A = 1.0, 2.3 
and 3.0, respectively; (b), K = 1.6, curves 1, 2 and 3 corre- 
spond to A = 1.0, 2.0 and 3.5, respectively. The chain length 
is TV = 4181. 

The Thouless exponent is proportional to the inverse 
of localization length i.e. 7 ~ l/£. If 7 is about order 
of 1/N for finite chain of length N, then the eigenstates 
are extended or critical. Otherwise the eigenstates are 
localized. In Fig. ||, we plot 7 as a function of eigenstate 
for different A. Fig. ||a is for K — 0.4. It tells us that for 
small A, all states are extended and critical (correspond 
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to small 7 in Fig. |]a). As A is increased, some eigen- 
states become localized. For A > A c all eigenstates are 
localized. Of course, A c depends on K @. The minimum 
Thouless exponent "f m in and maximum PR are calculated 
as a function of A so as to find the critical value of A c for 
a fixed K. In Fig. ||, we plot j m in an d PRmax versus A 
for K — 0.4, 0.6, K c , and 1.6. From these curves, we can 
easily obtain critical values A c « 2.3, 2.46 and 2.96 for 
K = 0.4, 0.6 and K c , respectively. 

In the case of K > K c , i.e., the ground state configui 
tion of atoms corresponds to Cantorus, both spectra a 
eigenstates are quite different from that case of K < , 
(see Fig. ||b). In this case, no critical value has be 



The variance of the wavepacket is 

N 

cr 2 (i)=£(n-n) 2 |^ n (t)| 2 . (8) 

n=l 

It can be calculated numerically by integrating the 
Schrodinger Eq. (Q) for a chain of length N with fixed 
boundaries tf>o — ipN+i = 0. In our calculations, the 
fourth-order Runge-Kutta method with time step 6t — 
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FIG. 3. The minimum Thouless exponent jmin (a) and 1 
maximum participation ratio PRmax (b) as functions of A 
electron in the FK chains with different parameter K. Cun 
1, 2, 3 and 4 in (a) correspond to K = 0.4, 0.6, K c and 1.6. 
Curves 1 and 2 in (b) correspond to K = 0.4 and 1.6. The 
results of 7 m i„ and PRmax are obtained for finite FK chains 
of length N = 4181 and N = 1597, respectively. 

To investigate quantum dynamical behaviors, the time 
evolution of a wavepacket in the system described by Eq. 
(Q) is calculated numerically. The wavepacket is localized 
initially at the centre of the chain. The time evolution is 
described by a time-dependent Schrodinger equation 



.dlp n 
I — ; — 

dt 



ip n+ i + ip n -i + Acos(x°)^« 



(7) 
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FIG. 4. The variance a 2 of a wavepacket for electrons in 
the FK chains with different values of K. (a), K = 0.4, curves 
1, 2, 3, 4 and 5 correspond to A = 1.0, 2.2, 2.3, 2.4, and 3.0, 
respectively, and A c ~ 2.3; (b), K = 1.6, curves 1, 2, 3 and 
4 correspond to A = 1.0, 2.0, 3.0 and 3.5, respectively. The 
chain length is N = 10946. 

In Fig. [|, we plot cr 2 (t) for several values of A. Figs. 
|]a and [|b correspond to K — 0.4 and 1.6, respectively. 
For K < K c , the time evolution of electron is also sim- 
ilar to that of the Harper model, that is a 2 ~ t 2 and t° 
for A < A c (A = 1.0,2.2) and A > A c (A = 2.4), respec- 
tively. The unbounded diffusion in the regime of A < A c is 
caused by the existence of extended and critical states as 
discussed above. However, at A = A c , the system displays 
anomalous diffusion behaviors, and the dynamical expo- 
nent depends on K, and A sa 1/4 for K — 0.4 (see Fig. 
|]a). It will be of great interest to connect the exponent 
of anomalous diffusion with the multifractal dimension of 
critical eigenstates. For K > K Cl the time behaviors of 
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wavepacket are similar to that of quasiperiodic systei 
i.e. the dynamical exponent A depends on A. 



Now we turn to level statistics of the system and 
relationship with the dynamical exponent. Geisel et. 

observed that for bounded uncountable set of levi 
it is possible to count the number of energy gaps lar; 
than s and to calculate the integrated level-spacing c 
tribution (ILSD) defined by Pi n t(s) = J p(si)dsi. 1 
derivative of ILSD p(s) = —dpi n t/ds gives probabil 
distribution of level spacings. In Fig. ||a we show IL 
at A = A c for K = 0.4 and 0.6. It can be seen that e< 
distribution consists of two parts: exponentially dec 
part for small s and power-law decay part for large s. 
is well known that the localization of eigenstates resr 
in the Poisson distribution in energy level spacing stai 
tics |2l[]. It is thus reasonable to attribute the expon 
tially decay to the localized states, and the power-law 1 
decay to the critical states, as is the case in the Har] 
model and the quasiperiodic model. Fig. |5|b shows IL 
for several values of A with K — 1.6. Each distributioi 
very similar to that of the Harper model at A c and tl 
of quasiperiodic models. Unfortunately, the power-] 
like behavior of ILDS is not very significant for largi 
due to limited bin sizes. 

The exponent (3 for level-spacing distribution defir 
by p{s) ~ s~P can be obtained by best fit power-law p 
of ILSD curves. We find that the relation 



A = 1 - j3, (9) 



for the Harper model at A c is also true for our model in 
the regime of K < K c , but not for the regime of K > K c . 

In summary, we have studied spectral and dynamical 
properties of an electron in incommensurate FK chains. 
The system shows rich phenomenon. For K < K c , i.e., 
the ground state configuration corresponds to invariance 
circle, there exists a critical value A c above which all 
eigenstates are localized. Below A c , extended, critical 
and localized states coexist. The critical value A c and 
dynamical exponent A at A c depend on K . The relation 
holds for power-law part of ILSD. On the other hand, 
for the case of K > K c , i.e., the ground state configura- 
tion corresponds to Cantorus, all electron eigenstates are 
critical and A depends on A. 
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FIG. 5. Integrated level-spacing distribution Pi n t{s) for dif- 
ferent values of K. (a), Curve 1 corresponds to K — 0.4 and 
A c = 2.3. Curve 2 corresponds to K — 0.6 and A c = 2.46; (b), 
K — 1.6, curves 1, 2, 3 and 4 correspond to A = 1.0, 2.0, 3.0 
and 3.5, respectively. The chain length is N = 2584. 
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